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Abstract 

Local quark-hadron duality violations in conventional applications of the 
operator product expansion are proposed to have their origin in the fact that 
the QCD vacuum or a hadronic state is not only characterized by nonvanishing 
expectation values of local, gauge invariant operators but also by finite corre- 
lation lengths of the corresponding gauge invariant n-point functions. Utiliz- 
ing high-resolution lattice information on these correlators a non-perturbative 
component of OPE coarse graining is proposed which, in principle, allows for a 
determination of the critical dimension where the break-down of the expansion 
sets in. 



A practical and very successful approach to relate hadronic parameters like 
masses, decay constants and widths to a few universal parameters of the vacuum 
of Quantum Chromodynamics |1| was introduced a long time ago The method 
of QCD sum rules (QSR's) relies on Wilson's operator product expansion (OPE) 
P[ at an external, euclidean momentum Q, analyticity of the considered QCD cur- 
rent correlator everywhere in the plane except for a cut along the positive real axis, 
and perturbative calculability of the operator coefficients, non-perturbative effects 
are introduced as power corrections via the nonvanishing vacuum averages of local, 
gauge invariant operators. If the average of an operator of mass dimension d is com- 
parable with A d , where A is a typical hadronic scale then modulo numerical factors 
power corrections are suppressed by (A/Q) d , and one may hope that the expansion 
approximates the correlator at least in an asymptotic sense ||. To make contact 
with hadronic properties one assumes quark-hadron duality, namely the property of 
the dispersive part of the correlator to be described in terms of measured hadronic 
cross sections. By appropriately adjusting the external momentum QSR's relate the 
properties of the lowest resonance to the non-perturbative condensates. There is not 
much doubt that averaged spectral functions relevant in QSR's are dual to the OPE 
0| . In practice it is fair to say, however, that there are channels where the sum rule 
method seems to be jeopardized by abnormally large non-perturbative corrections 
(scalar, pseudoscalar) . 

However, it is painfully obvious that the conventional, practical OPE, which is 
typically truncated at dimension d = 6 or d = 8 does not account for pointwise or 
local duality (for a summary see the review articles ||, |5[|). Usually the asymptotic 
nature of the OPE is blamed for this. Indeed, perturbatively calculated Wilson coef- 
ficients being of the form (a s (Qo)/(«s(Q)) 7 x (Q 2 )~ d , where 7 denotes an (effective) 
anomalous dimension and Qo is a normalization scale, do not allow for the observed 
"wiggling" of the spectral function obtained by taking the imaginary parts of a 
term-by-term continuation to time-like momenta. Since the expansion is truncated 
terms, which behave like exp [—const J Q 2 / A 2 ] , can not be captured ||. On the other 
hand, such terms contribute asymptotically unsuppressed oscillations to the spectral 
function which contradicts experiment and is not allowed by asymptotic freedom. 
Investigations concerning the origin of local duality violation have been performed 
within various models in the literature H |7|, |8|. Impressively, it was shown in |5], I7J 
that upon appealing to a dispersion relation a model spectral function with equally 
spaced narrow resonances yields an asymptotic expansion in the euclidean which 
resembles a conventional OPE with factorially growing coefficients. However, such a 
model ignores the drastic broadening of resonances and the decrease of their heights 
with increasing energy. 

Local duality violations may have a considerable impact on the calculation of 
life-time differences and inclusive/exclusive decay widths of S-mesons where OPE 
in conjunction with the heavy quark expansion is applied at the large time-like 
momentum p ~ rrib || |10[ . The precise theoretical determination of these quantities 
is of particularily acute relevance since it would lead to stringent constraints on the 
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parameters of the underlying model (e.g. CKM matrix elements in the Standard 
Model (SM) and/or couplings, masses, and CP violating phases of SM extensions). 
This, in turn, is of paramount importance for the potential detection of New Physics 



It is therefore essential to improve our understanding of local duality and to 
estimate the scales which govern it. 

The purpose of this letter is to point out and exploit a genuine source for OPE 
modification which has only been vaguely addressed so far ||. The observation is 
that the conventional OPE of a current- current correlator, which can be proven to 
exist perturbatively, does not account for the diminishing relevance of local opera- 
tors composed of the fundamental fields in an asymptotically free theory like QCD 



when decreasing the spacetime resolution [T^|. This resolution, however, is deter- 
mined by the external (euclidean) momentum squared Q 2 = —q 2 . The physical 
picture behind this can be examplified in the case of a correlator of light quark 
currents as follows D: An external current creates a quark-antiquark pair with 
total momentum q. Due to asymptotic freedom the propagation of these degrees 
of freedom is dominated by perturbation theory at high enough Q = y/Q 2 and 
comparable momenta, non-perturbative effects become important if there are soft 
emissions in the process, which are effectively described by nonvanishing vacuum 
expectation values (VEV's) of gauge invariant, local operators. But what does lo- 
cal really mean? We believe that the answer to this question is tightly connected 
to the issue of pointwise duality. During the time interval between emission and 
reabsorption of the constituent degrees of freedom of a hadron the vacuum, which 
admits these propagations, is being resolved precisely by the external momentum 
Q. So if the OPE is an expansion involving local operators one also ought to say 
how resolved the corresponding locality is. The renormalization group teaches that 
information on the VEV of some local operator at a given resolution derives from 
local and non-local information at higher resolution. In perturbation theory, where 
interactions do not generate effective fields and parameters out of their fundamental 
counterparts, this is done by normalizing the corresponding operator at a definite 
momentum scale Qo and relating this to the momentum Q of the process by means 
of the corresponding set of renormalization group equations. Since in QSR's the 
VEV of a local operator is supposed to exclusively carry non-perturbative informa- 
tion one must apply an exclusively non-perturbative procedure of coarse graining to 
it. The key to non-perturbative coarse graining for the lowest dimensional operators 
is information on the corresponding gauge invariant 2-point functions [|I3[] which is 
available from lattice calculations [fPjl . For higher dimensional operators one so far 
has to appropriately impose vacuum saturation to express their correlators in terms 
of the 2-point functions. 

The idea that due to the existence of finite correlation lengths in gauge invari- 
ant 2-point functions non-local, non-perturbative effects should be included in the 
OPE is not new. We just mention here that to obtain realistic hadronic wave func- 
tions along the lines of ref. [0 non-local effects were systematically included via an 
expansion in covariant derivatives in ref. [JT7]. On the other hand, based on pertur- 
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bative arguments it was shown in refs. |18| that a light-cone expansion of the product 



of gauge invariant currents can be written in terms of non-local string operators. 
Working in position space, perturbative evolution equations were derived for these 
operators in In this paper we introduce non-perturbative, non-locality within a 
phenomenological generalization of the standard Wilson OPE. Although we concen- 
trate on vacuum correlators it is clear that the procedure is extendable to hadronic 
state averages. Applications of the concept outlined above to OPE's and spectral 



functions of various light quark channels are reported in ref. ||20|| . 

We start out by expanding the correlator into a conventional OPE at a large 
external, euclidean momentum Q, where we expect the description of the dynamics 
in terms of the continuum action and local operators made of fundamental fields to be 
sufficiently accurate. To proceed to lower external momenta Wilson coefficients are 
evolved perturbatively by means of the running coupling a s (Q 2 ) and the anomalous 
dimensions of the operators. The evolution of an operator average can be obtained 
using information on the non-perturbative part of the corresponding gauge invariant 
correlator in euclidean position space 

(F(0))y = (Fi(0) ■ • -F B (0))J - l/NY. <*i(0) ■ ■■F n {x n )Y* . (1) 
Thereby, appropriate parallel transporters 



S(0,x) = Pexp 



ig / dz n A n . (2) 
J o 



in the non-local expression are understood to make the object gauge invariant. The 
sum runs over all possible (piecewise straight trajectories of parallel transport. 
A normalization 1/M depending on the mass dimension d and the numbers of fields 
transforming under the fundamental and adjoint representation is implied to re- 
duce the correlation function to the "condensate" in the limit x%, ■ ■ ■ , x n — > 0. As 
examples we give the following dimension 5 and 6 operator VEV's 

(z) (^(0)F^(0K^0))J, iii) (trF^(0)F, K (0)F KM (0))J p . (3) 

If an arrow pointing from X{ towards Xj stands for the parallel transport S(xi,Xj) 
then the relevant contribution to the correlators correponding to the local operators 
of eq. (||) can be depicted as in Fig. 1. Note that points with fields transforming 
under the fundamental (adjoint) representation are connected to one (two) lines of 
parallel transport. 

To coarse grain from resolution Q to resolution Q — dQ we propose to average the 
non-perturbative part of a correlator corresponding to a connected diagrams over a 
(euclidean) ball of radius dRQ with 



1 11/1 \ dQ 
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Figure 1: Diagrammatic representation of the gauge invariant correlators corre- 
sponding to the local operators of (i) and (ii) in eq. (H). 



In the short-hand notation of eq. (|l|) this is written as 



1 

(V{dRQ)f n - 1 ) J\ Xl \,..,\ Xn \<dR Q 



(m) n Q P - d Q = n„ Jn uT^l / _ d 4 x 1 ---d*x n (F 1 (0)...F n (x n )) n o p , (5) 

where 




V(dR Q ) = \^{dR Q f = ^ ( ^ ) . (6) 

As for correlators described by disconnected diagrams the application of the fac- 
torization hypothesis (vacuum saturation) is discussed in |2D| for the example of 
4-quark operators composed of color singlet or octet currents. 

Let us now concentrate on 2-point functions as they are relevant for dimension 
3 and 4 quark and gluon operators, respectively. To our knowledge only the gauge 
invariant bilocal quark correlators (vector and scalar) |TJJ and the gluonic field 
strength correlator |L3| have been measured on the lattice. The results suggest that 
there exists an additive decomposition into a perturbative, power-like in \x\, and a 
non-perturbative, exponential in \x\ piece fT^Q We are interested in the latter 

(F 1 (0)F 2 (x))% = A(Q)e*p(-\x\/\). (7) 

Since A is expected to have direct physical meaning - see for example ref . [jnj] , where 
in the case of the field strength correlator A -1 was interpreted as the energy difference 
A between a heavy hybrid (an adjoint quark with glue) and the mass of the heavy 
quark - it can not depend on the resolution Q. So coarse graining may only affect the 
pre-exponential factor A(Q). To derive an evolution equation for A(Q) we combine 
eqs. (|) and (0): 

(v(dQ/Q 2 )) 1 [ rfn d 4 x A(Q)e'^ x = A(Q - dQ)e~ / x = A(Q - dQ) . (8) 



1 There have been a number of attempts to calculate gauge invariant quark and gluon correlators 
in the instanton vacuum p5[. In this work, however, we do restrict ourselves to direct lattice results. 
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Note that in the limit A — > oo this leaves A invariant. In this limit we recover 
the conventional treatment of non-perturbative corrections in the framework of the 
OPE. Expanding the l.h.s. and r.h.s. of eq. (Q) in ^ and comparing coefficients of 
the linear terms, we obtain 



(9) 



Comparison of zero order terms yields an identity as it should. Evolution equations 
like eq. (H) could be derived for (n > 2)-point functions if the explicit form of the 
non-perturbative part could be extracted from the lattice. The solution of (pF) is 



A(Q) = A(Q ) exp 













Qo)_ 



(10) 



Disregarding powers of logarithms in the Wilson coefficients, which come from the 
anomalous dimensions of the corresponding operators and combined with the run- 
ning coupling a s (Q 2 ) possess exponents which are numerically small for the naively 
vacuum saturated 4-quark operators in the OPE's of light-quark correlators [Q, 
eq. (|ToD implies the following, generic form for the non-perturbative corrections of 
dimension 4 and 6 



(Q 2 ) 2 



exp 









5A 4 




Qo)_ 



A 6 (Qo) 

{Q 2 f 



exp 




Here the term 'naive' refers to vacuum saturation at the level of local operators and 
subsequent derealization of the scalar 2-point functions (see also ref. P0fl ). Eq. ( pT)| ) 
implies that conventional, dimension 4 and 6 power corrections are drastically altered 
if Q is of the order of A -1 or less. Despite the fact that the expression for dimension 
6 in eq. (|ll|) heavily relies on the vacuum saturation hypothesis it may feature a 
general trend of decreasing (effective) correlation lengths in higher mass dimensions 
(ultimately, only the lattice can decide on the truth of this assertion). 

Let us extrapolate the behavior of eq. (|Tl"D to non-perturbative corrections of 
arbitrary dimension d. If the running with Q essentially can be expressed as 



A d (Q ) 




(12) 



where A^ denotes an effective correlation length, then there is a single maximum at 
momentum Q = (dXd)^ 1 with value 



P d = (A d X d exp[-l] d) d exp 



a; 



Qc 



(13) 



where = {A d {Qo)Y l ^ d \ So a non-perturbative correction at fixed d vanishes 
asymptotically for both, Q — > oo and Q — > 0, which is in stark contrast to the 
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conventional picture. This, however, does by no means imply that we approach 
asymptotic freedom at low momenta since higher and higher mass dimensions may 
contribute their maxima at lower and lower Q. One possibility for this to happen is 
the following: Let Ad fall off as A^ = \/S- l ~ £ \ (1 > e > 0), then the positions of the 
maxima still decrease as Q = \~ l d~ e . The values of these maxima behave as 

P d = (A d Acf exp[-l + ^-]^ . (14) 

In particular, assuming that = Aqcd as it is suggested by conventional OPE's, 
the critical dimension d c from which on maxima explode is given implicitely as a 
function of A, e, Aqcd, and Q by 

A QCD Ad £ exp[-l + ^] = 1 • (15) 

Truncating the OPE at some d t < d c would mean that physics below the corre- 
sponding Q is not described. A truncation at d t > d c does not make sense since due 
to the rapid increase of the contributions at d > d c we are surely at orders where 
the expansion does not approximate (asymptotic expansion). Let us give some nu- 
merical examples. For e = 0.5, Aqcd = 0.4 GeV we calculate (the rounded) d c as a 
function of A 

A = 3 GeV" 1 -> 4 = 4; A = 2 GeV" 1 -> 4 = 10; 

A = lGeV" 1 -> 4 = 40; A = 0.5 GeV" 1 -> 4 = 158. (16) 

Lattice measurements with Np = 4 staggered fermions of mass a ■ m q = 0.01 and 
lattice resolution Qq = a" 1 ~ 2 GeV suggest that fermionic and gluonic correlation 
lengths corresponding to scalar correlation functions are ~ 3.2 GeV -1 and ~ 1.7 
GeV" 1 , respectively 0]. So for dimension 4 and assuming (naive) vacuum satu- 



ration also for dimension 6 this does not pose a problem for the conventional sum 
rule analysis of light quark correlators (at Qo = a" 1 ~ 2 GeV Ag q (Q ) as well as 
Af 2 (Qo) seem to be compatible with their phenomenological values at \i ~ 1 GeV 
lUll). However, (naive) vacuum saturation and/or the exclusion of higher mass di- 
mensions do not lead to a realistic behavior of the spectral functions which is shown 
in [20]. Still, going only to dimension 6, one obtains a much more realistic pic- 
ture if the (effective) correlation length at this dimension is drastically decreased 
(by about 10 times!). This can be motivated by the small correlation length of the 
2-point vector function on the one hand (naive vacuum saturation not applicable). 
On the other hand this short correlation length may cumulatively simulate higher 
mass dimensions. The tendency that dimension 6 contributions are overestimated 
in the conventional OPE of 10 for distances > 0.3 fm can be extracted from the 



work of p2[ . There, the euclidean position space V ± A correlators were calculated 



in a random instanton liquid model and compared to the ALEPH data [Elf and the 
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conventional OPE. Note, however, that fitting experimental data to a conventional 
OPE does violate local duality. We investigate the issue in p3 . 



In the modern literature there is supporting evidence in favor of large mass scales 
(as compared to Aqcd ~ 0.2 — 0.4 GeV) characterizing the QCD vacuum: Based on 
lattice data for the instanton size distribution it was stated in ref. Shuryak that the 
average size of instantons is much smaller (~ 1/3 fm) than their mean separation 
(~ 1 fm). Using maximal Abelian projection on the lattice to identify magnetic 
monopoles, the size of these defects, r ~ 0.05 fm, was determined for an SU(2) pure 
gauge theory (/3 ~ 2.5) in ref. |26]. To calibrate the lattice the SU(3) value of the 
string tension a = 440 MeV was used. It was found in p7) that a phenomenologically 
motivated tachyonic gluon mass can take values well above the 1 GeV bound. 

To summarize, we proposed to use non-local, non-perturbative information on 
the QCD vacuum (extendable to hadronic states) for a non-perturbative procedure 
to coarse grain the euclidean VEV's (or hadronic matrix elements) of local operators 
composed of fundamental fields down to resolutions where their relevance must be 
questioned. Implications for the low-energy behavior of the OPE in the euclidean 
were discussed. In particular, a possibility for the determination of the critical 
operator dimension, where the breakdown of the OPE sets in, was outlined. 
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